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ABSTRACT 

From a bibliographic analysis of the major published research studies about the index values used to evaluate 
water quality, the general characteristics of the so-called Biological Monitoring Working Party (BMWP) index are 
described; the degree of complexity of water quality assessments is measured and the integration of new technologies is 
analyzed through the use of Geographic Information Systems (GIS). The BMWP index is studied, as well as its 
applications from 1978 to 2017, specifying its temporal and spatial distribution, its current trends and complementary 
values. The analysis is based on the application of the Graph Theory, particularly on the Social Networks Analysis (SNA) 
from a holistic approach. 
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INTRODUCTION 

The quality of water resources is altered by anthropogenic activities, population growth, and natural 
phenomena, among other factors, which stress ecosystems and risk water availability. The approach to assess it has 
been developed on three fundamental pillars: physical, chemical and biological. For this purpose, the use of 
biological metrics stands out by being less expensive, less complex and depending on the bio indicator used they 
provide spatial and temporal information. (Barbour, Gerritsen, Snyder, & Stribling, 1999; Kumar, Colton, Springer, 
& Trama, 2013; Mangadze, Bere, & Mwedzi, 2016; Soldner et al., 2004; Yazdian, Jaafarzadeh, & Zahraie, 2014) 
Moreover the Biotic indices to assess water quality in freshwater ecosystems have been applied since last century 
(Armitage, Moss, Wright, & Furse, 1983; Hilsenhoff, 1977; Kolkwitz & Marsson, 1909; Woodiwiss, 1964) and its 
development with theoretical mathematical support (Margalef, 1951; Shannon & Weaver, 1998; Simpson, 1949) 
became popular in the 50’s. 

The Biological Monitoring Working Party index (BMWP) takes into consideration the presence or absence 
of aquatic macro invertebrates to evaluate water quality. It was developed based on the sensitivity of these family 
groups to different degrees of organic pollution in rivers, without taking into account the abundance of each taxon 
(Hawkes, 1998). It grew to be used as a worldwide reference (Basset et ah, 2012; Damanik-Ambarita et ah, 2016; 
Forio et ah, 2017; Munoz-Riveaux et ah, 2003; Perez, 2003; Romero et ah, 2017; Rosas-Acevedo et ah, 2014), even 
though the original theory hadn’t been published. 
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In Spain, Alba-Tercedor and Sanchez-Ortega (1988) modify the tolerance values of each macro-invertebrate 
family, adapt the BMWP references to local conditions (BMWP’) and define five types and colors associated with the 
index value, to describe water quality, and in order to develop biological maps. On the other hand, the BMWP-Cu index 
(Munoz-Riveaux et al., 2003), is a direct consequence of the BMWP’. Both methodologies describe how to map water 
quality, but do not show any representation. Additionally, Perez (2003) suggests that the methodological theory for 
regionalization of water quality should be tailored to each region. On the other hand, Forero, Longo, Ramirez, Jairo, and 
Chalar (2014) in Colombia, develop the Multi-Metric Ecological Quality Index (ICE rn _ mae ), as a biological and physic- 
chemical method for the evaluation of water quality. It exploits the physical and chemical parameters, environmental 
gradients, abundance, and macro-invertebrate genus identification. Furthermore, they compare it with the BMWP-Col 
index (Perez, 2003) and prove that both indicators provide similar information on the variability of water quality. 

Kohlmann, Russo, Itzep, and Solis (2010) employ a cooperative methodology to assess water quality: using the 
BMWP’-Cr reference in two rivers of rural Costa Rican communities. Moreover, they integrate the participation of the 
scholar community through a survey, before and after the study of water quality. They conclude that the participants 
changed their minds about how they perceived their water-related environmental problems. On the other hand, the impact 
of anthropogenic activities on water quality in Iran is studied through the use of Geographic Information Systems (GIS) 
and comparative analysis between physicochemical and biological indexes -one of them the BMWP- and show that the 
biological indexes provide a better approach to the health level of ecosystems. (Sharifinia, Mahmoudifard, Namin, 
Ramezanpour, & Yap, 2016) 

In this way, a complete water quality assessment is governed by the monitoring of three major components: 
hydrological, physicochemical and biological. The latter, should target the response of species or communities to changes 
in their environment, and determine what the factors and activities that interfere directly or indirectly on water quality. 
(Chapman & Organization, 1996; Hellawell, 2012) Moreover observed that, over time, water quality studies using the 
BMWP index, have also been integrating: new technologies, the participation of social actors or anthropic activity, more 
complex relationships between physicochemical, biological and environmental variables, the index itself as an aid to other 
metrics and their adaptation in lentic ecosystems(Castellanos Romero et al., 2017). 

This article attempts to give a holistic approach to water quality studies, based on the BMWP index methodology, 
through measuring the degree of complexity of water quality studies, the application of the index in different ecosystems, 
the incorporation of new technologies (GIS) and the use of other metrics in water quality assessments. 

MATERIALS AND METHODS 

The bibliographical sources consulted present aquatic macro invertebrates as bio-indicators for the study of water 
quality, besides the application and adaptation of the BMWP index in different parts of the world, from 1978 to 2017. 
Scopus and Web of Science were used to identify the most relevant studies in specialized journals, in both English and 
Spanish. The research began with the keywords “BMWP index’’ and “water quality”, from which 177 and 119 
bibliographic sources respectively were obtained, as well as 25 more results from ‘Google Scholar’. In the next phase, 
databases were integrated to avoid duplicate sources, and 45 records were discarded, given that their results were not 
relevant to the study. Among the 165 remaining studies, a random sample of about 62% (98 articles) was chosen. 
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Table 1, describes the set of variables that are used for each particular aim. The data matrix is shown in Appendix 
1. The degree of complexity was measured with the number of parameters used in each article consulted and is described in 
Table 2. To determine if the number of parameters used over time increased or decreased, the Pearson correlation 
coefficient was calculated to explain the relationship between the year and the number of parameters (IBM SPSS Statistics 
ver. 20). Google Earth 7.3.0 was used for geo-referencing rivers, streams, reservoirs, lakes, and other water bodies, in order 
to describe the spatial and temporal distribution of the use and adaptations of the BMWP index; the graphic representation 
was made with ArcMap 10.3. On the other hand, «cartography and GIS» includes four variables that are described in 
detail in Table 3. 


Table 1: Objectives Description and Their Relation with Each Variable Defined 


Objective 

Variable 

Description 

Measure the number of parameters used in 
water quality studies. 

Sources and parameters (see, 
figure 2) 

Author and physical, chemical and 
biological parameters 

Explain the relationship between the year and 
the increase in the number of parameters used 

Year and parameters 

Source publishing year and 
parameter numbers measured in 
the bibliographic sources 

Identify the BMWP adaptation / application 
trends at an Ecosystem level 

Source, ecosystem and year 

Author, lentic or lotic ecosystem 
and source publishing year 

Describe the spatial and temporal distribution 
of the BMWP application / adaptation 

Country and year 

Country and source publishing 
year 

Detail the cartography and GIS used in water 
quality 

Source, year and, 
«Cartography & GIS» 

Author, source publishing year 
and see Table 2 

Identify the indexes or metrics used in the 
assessment of water 

Several Indices 

Biological, chemical and 
multimeter 


Since it is possible to obtain quantitative data or significant statistics from the collection of the data described 
here, in this review we used the Social Network Analysis (ARS), which includes a set of qualitative measures to describe 
the structure and development of networks, which allowed describing the general relationships between the variables in 
Table 1, in addition, a graph was developed with Gephi 0.9.1 (Bastian, Heymann, & Jacomy, 2009). The graph analyzed is 
not directed (simple), with nodes (vertices); associated with the author, ecosystem, country, year, parameters, 
«cartography and GIS» and indexes. The edges represent the relationships: Source/parameters, source/ecosystem, 
source/«cartography and GIS», country, year and indexes (Wasserman & Faust, 1994). In correlation, the modularity 
algorithm was applied (Blondel, Guillaume, Lambiotte, & Lefebvre, 2008) to describe the similarity between communities. 

RESULTS AND DISCUSSIONS 

Figure 1, describes the variables behavior in a global way. The revised studies are arranged as central nodes 
(author or source), labeled with a number, sorted down from the highest node degree (the one with the most correlations). 
The index has been applied in lotic ecosystems. 2014 is the year with the most published sources. The indexes BMWP, 
ASPT, Shannon-Weaver, and Margalef are the most employed. On the other hand, 61.2% of the consulted research studies 
use cartography or satellite images for the study area description. In addition, the most used parameters were: Dissolved 
Oxygen, Temperature, pH, Electrical Conductivity. 

ARE WATER QUALITY STUDIES COMPLEX? 

Among the 98 reviewed research studies, only 42 parameters were considered for the source-parameter 
correlation, but others such as latitude, altitude, slope, substrate, and riparian vegetation were not considered. However, it 
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should be noted that some of these parameters are analyzed within the Cartography or GIS section. The ‘parameter’ 

variable showed values between 0 and 23 as minimum and maximum. The average value of measured parameters was 

7.3+5.58 per source consulted. The quartiles values were: 2.75 (First), 7.5 (Second) and 11 (Third). 

On the 19% of the consulted sources did not report some parameter, since the selected sample did not discriminate 
between studies that diagnose water quality and methodologies that explain the adaptation of the index. It is wise to 
recognize there are methodologies that include case studies where the proposed index is validated. On the other hand, 4% 
of the reviewed sources measure more than 18 parameters, which indicate that water quality assessment is more rigorous in 
some countries. According to the quartiles, the results of water quality studies are classified and shown in Table 2. 

Dissolved oxygen, pH and temperature are the most measured parameters, and the match the top three records 
sent to the Water Quality Monitoring Data Repository of EPA. gov (STORET), and they’re used to determine the 
minimum set of parameters that must be measured, to consider a reliable study of water quality, while for particular or 
contextualized studies other parameters should be included. For example, total and fecal coliforms, arsenic, mercury, 
boron, and, nitric acid, among others. 
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Figure 1: Graph of the Data Matrix, Where Variables Grouped By Colors 
Table 2: Classification of Quality Studies According to the Measured Parameters 


Parameters 

Study 

Frequency 

Percent 

0-3 

Diagnostic 

27 

27.6 

4-8 

Basic 

31 

31.6 

9-11 

Half complex 

18 

18.4 

>= 12 

Complex 

22 

22.4 
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Year 


Impact Factor (JCC): 4.6869 


Figure 2: Scatter Diagram of the Year/Number of Parameters 
Correlation (n=98 y p=0.05) 
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To explain water quality study complexity in recent years, let us analyze the scatter diagram in Figure 2. It shows 
the correlation between the research publishing year and the number of parameters measured. The Pearson correlation 
coefficient was R : = 0.081, proving that there is no positive linear correlation between the year of the source and the 
increase in the number of parameters measured. This indicates that the studies carried out to evaluate water quality are 
highly complex and, essentially, contextualized. Basically, it is because North America and the Western Asia countries are 
still developing basic methodologies related to the BMWP index. In Figure 2, the points enclosed within the rectangle 
belong to those water quality studies carried out in the last 11 years, which only applied the BMWP index or measure less 
than 11 parameters. 

SPATIAL DISTRIBUTION AND TRENDS 


Regarding the origin of the sources, the published research studies from Brazil, Colombia, Spain, and the United 
Kingdom constitute more than a third part of the consulted resources. The rest of the papers are distributed in 33 countries 
from around the world. Figure 3, shows the 37 Source/Country communities, in which the quantity and origin of the 
research articles reviewed are shown. The community size is determined by the number of nodes; the degree of the nodes 
represents the number of papers reviewed and the label shows the country of origin for the research. 

For a more detailed analysis of the spatial and temporal distribution of the BMWP index application, the study 
was filtered by decade. From 1978 to 1990, the BMWP index was applied in the United Kingdom and Spain. Its use 
increases from 1991 to 2000, mainly in Europe, in countries like Portugal, Poland, Greece, and Slovakia. In addition, in the 
same period, its application reaches Latin American countries such as Argentina, Brazil, and Colombia. From to 2001 to 
2010, it starts getting used in Africa (Egypt and Mozambique), Asia (Iran, Malaysia, and Thailand) and other countries in 
Latin America (Cuba, Costa Rica and Venezuela, among others). From 2001 to 2010, it began its application in Africa 
(Egypt and Mozambique), Asia (Iran, Malaysia, and Thailand) and other countries in Latin America (Cuba, Costa Rica, 
Venezuela and among others). From 2011 to 2017, Mexico, Russia and other countries in Asia (India and Vietnam) share 
its use. This period is the one with the most water quality studies reported (44). A general description is shown in Figure 4, 
elaborated with data obtained from the geo referencing of the reviewed papers. 



Figure 3: Community Distribution by Country/Source. The Graph on the Right Shows the 
Node Number in Each Community. The Red Circles Represent the Largest Communities 


Note that in Europe there are a large number of water quality studies, which are based on the absence or presence 
of macro invertebrates and the BMWP index, or any of its variations, followed by America, Asia, and Africa. Since its 
initial theory, the BMWP index has been applied to lotic ecosystems, but since 2013 its application in lentic ecosystems 
has increased. 
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CARTOGRAPHY AND GEOGRAPHICAL INFORMATION SYSTEMS 

Among the 98 reviewed papers, only five research studies were conducted in lentic ecosystems areas, 
corresponding to countries such as Colombia, Ethiopia, Guatemala, India, and Mexico. This study shows that only three 
countries (Colombia, Mexico, and Ethiopia) study the water quality in both ecosystems. In this way, the popularity, 
adaptability and reliability of the BMWP index is manifested due to the fact that after 38 years since the developing the 
BMWP index, 44.89% of the studies have been carried out in the last six years, and that in the last four years it has been 
applied to other ecosystems. 
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Coordinate System: GCS ETRS 1989 
Datum: ETRS 1989 
Units: Degree 


Figure 4: Spatial and Temporal Distribution of the BMWP Index 
Application or Adaptation (1978 - 2017). 


Four variables were classified to evaluate the «Cartography and GIS». Their description and nomenclature are 
shown in Table 3. 


Table 3: Variables that Describe the Cartography and GIS Use 


Variable 

Description 

Study area (SA) 

Satellite image or map that only describes the study area. 

Water Quality (WQ) 

Cartography that uses intervals and water quality classes, according to 
(Alba-Tercedor and Sanchez-Ortega (1988)). 

Landsat (LT) 

Satellite image that describes the study area and other environmental gradients (GIS) 

LandsatWQ (LTWQ) 

Satellite image used to relate environmental gradients with water quality (GIS). 


Among the consulted sources two large groups can be spotted, among which 24 analyze the study area 
descriptively (zero degrees node) and 60 describe the study area with maps or satellite images. Eight use satellite images or 
maps to study the area and other environmental variables. The rest of them are defined in Figure 5: the red circle encloses 
the revised sources that use water quality cartography as described by Alba-Tercedor and Sanchez-Ortega (1988). 

Of the sources analyzed, only 4.9% use water quality cartography, maps, and satellite images to describe the study 
area. And only two integrate the four variables described in Table 2. The fact that only five of the 98 sources reviewed use 
regionalization by colors described since 1988 by Alba-Tercedor and Sanchez-Ortega (1988) and later by Perez (2003) and 
Munoz-Riveaux et al. (2003); shows the absence of an explicit methodology for the elaboration of water quality maps 
associated with BMWP index. Although in more recent Works (Romero et al., 2017) the use of GIS for the regionalization 
of water quality is observed, without describing the methodology used for the elaboration of maps of water quality, 
necessary in the spatial and temporal description of the conditions of the body of water. 
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Figure 5: General Correlation Cartography/Sources 


Parameters such as latitude, altitude, slope, depth, chlorophyll, substrate, riparian vegetation, temperature, total 
suspended solids and nitrates can be calculated through satellite images, by the use of multispectral sensors (Bierman, 
Lewis, Ostendorf, & Tanner, 2011; Gomez Diaz & Dalence Martinic, 2014; Norat Ramirez, Mattei, Barreto, & Luhring 
Gonzalez, 2002; Rivera Vazquez, Vargas Perez, Terrazas Dominguez, & Gavi Reyes, 2002; Trangmar, Yost, & Uehara, 
1986). For example, depth is related to bathymetry, riparian vegetation is related to vegetation or type of land use. 
Furthermore, nowadays there are digital elevation models that include: latitude, altitude, and slope. 

OTHER METRICS OR INDEXES 

The 277 nodes shown in Figure 6 are distributed in 142 indexes (physical, chemical, biological and other metrics), 
98 revised research study papers and 37 countries around the world. The degree of the nodes associated with BMWP, 
BMWP’, ASPT, Shannon-Weaver, Margalef, EPT, WQI and Simpson indexes are the most cardinal in the graph studied, 
which shows that they are the most widely used. This implies the transition from the qualitative to a quantitative sampling 
type to describe the ecosystems health state. Besides, the application of the WQI index reflects on the presence of 
physicochemical parameters. 



Figure 6: Communities Composed of Nodes: Indexes, Sources, 
and Countries, With a Modularity of 0.515 


The 10 communities (modularity class) described in Appendix 2, reveals the similarities among the index 
application in different countries, and their size is determined by the number of nodes that confirms it. Countries like 
Egypt, Nigeria, Venezuela, Iran, Turkey, Brazil, and the United Kingdom, are part of the larger community (61 nodes), 
whilst Cuba, Panama and Guatemala are in the smallest, highlighting that Cuba uses more indexes than Panama and 
Guatemala for their water quality studies. 
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From the analysis of papers consulted, associated with evidence, the importance of the contextualization of water 
quality studies is demonstrated. Therefore, it is necessary to add new taxa or to adapt the tolerance values of the macro¬ 
invertebrate families of each region and ecosystem, for the BMWP index methodology application. This reinforces the idea 
of the BMWP adaptability and validation. A sample of the foregoing appears in the bibliography as Y-BMWP, BMWP- 
THAI, BMWP-Viet, BMWP/Atitlan, BMWP'/PAN, BMWP-CETEC, BMWP (PL) y BMWP-Coyuca. 

On the other hand, index application (node degree less than 3) is relative to the complexity of the particular 
studies and the type of correlation to be analyzed. Moreover, the increase in emerging pollutants makes it necessary to 
resort to other metrics in order to have a better perspective the health of the ecosystems. The large variety of indexes or 
metrics used in the evaluation of water quality means that there is no index that could be adopted in general, or that could 
work correctly in different ecosystems. 

It is known that water quality studies that use bio indicators, as a rule, are governed by variables such as study 
area, physical-chemical parameters, biotic indexes and correlations with other metrics. According to that, water quality 
studies can be described in relation to the variables observed, such as Parameters: physical, chemical, microbiological and 
environmental. They allow the description of particular elements and conditions. Indexes: biological, physicochemical and 
multimetric. They describe the ecosystem dynamics qualitatively and quantitatively and use one or several parameters for 
its calculation. GIS help describes the environment of water bodies to integrate other environmental variables, in addition 
to spatial and temporal analysis. 

On the other hand, some of the analyzed sources integrate aspects such as Anthropogenic activities, the activities 
that surround the water body or for which the water resource is used (Blijswijk, Coimbra, & Gra£a, 2004; Duran, 2006; 
Garcia-Criado, Tome, Vega, & Antolin, 1999; Gray & Delaney, 2008; Guevara Mora, 2011; Holguin-Gonzalez, Everaert, 
Boets, Galvis, & Goethals, 2013; Mesgaran Karimi, Azari Takami, Khara, & Abbaspour, 2016; Stein, Springer, & 
Kohlmann, 2008; Yoshida & Rolla, 2012). The inclusion of public and educational institutions (Rico, Rallo, Sevillano, & 
Arretxe, 1992). Study of invasive species (Arbaciauskas, Visinskiene, & Smilgeviciene, 2011; MacNeil & Briffa, 2009) or 
with the presence or absence of other ones (Bedford, 2009). Ecological evaluations (Sporka, Pastuchova, Hamerllk, 
Dobiasova, & Beracko, 2009), correlations with changes in the landscape (Yoshida & Rolla, 2012) and land use (Munyika, 
Kongo, & Kimwaga, 2014). 

The General Systems Theory (Von Bertalanffy, 1969) raises that the phenomena must be studied as a whole, that 
is, all the parts, the associated processes, the organization and the order that links them. This theory proposes the well- 
know holistic approach for the study of systems, taking into account the elements that intervene and correlations 
established between them. Hence, the ecosystem seen as a system allows interrelating the variables and aspects described 
in this review, to give a holistic approach to water quality studies. Figure 7 describes the general correlations in water 
quality assessment which apply the BMWP index (according to the grayscale, they go from basic to complex, also 
considering the increase in the number of variables studied). 

CONCLUSIONS 

The aim and scope of the research works that assess water quality through aquatic macro-invertebrates, which 
apply or adapt the BMWP index are heterogeneous and multidimensional. The parameters mostly used for this analysis 
(dissolved oxygen, pH and temperature) are in order of importance, depending on the direct impact on the communities or 
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species and the information that it’s contributes of the health of the ecosystem of immediate way. There is no correlation 
between the increase of the parameters studied and the year of study. The spatial and temporal distribution of the studies 
carried out in different countries that apply or adapt the BMWP index, from the theory to the initial application, let us 
certify that the United Kingdom, Spain, Brazil, and Colombia are the nations that contribute most of the scientific work in 
the field of water quality evaluation. In addition, this distribution shows a clear tendency of typical taxonomic groups to 
adapt to lentic ecosystems in recent dates. 



Figure 7: Water Quality Studies Viewed as a System 


There is not the methodology that describes in a general way the making of maps associated with water quality 
studies that use the types and BMWP index value as reference elements. It is evident that there is a tendency to use 
technology in these studies. For example, some parameters that used to be measured with traditional methods before, have 
been proven to be described through the use of satellite images or sensors nowadays. Among the 142 identified indexes, 
the most used in the specialized literature are the following: BMWP, ASPT, BMWP', Shannon-Weaver, Margalef, and 
Simpson. This corroborates that the most used sampling methods in the present for water quality assessment are varied and 
non-directional, that means they are both qualitative (BMWP, ASPT, BMWP') and, quantitative (Shannon-Weaver, 
Margalef, and Simpson) and the taxonomic identification at different levels (family, gender or species) is used in both of 
them. 
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